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The polarization dependence of nonequilibrium transitions in a multistable cavity-polariton sys- 
tem is studied under a nanosecond long resonant optical excitation at the normal and magic an- 
gle incidences with various polarizations of the pump beam. The temporal correlations between 
the frequency, intensity, and optical polarization of the intra-cavity field, which all undergo sharp 
threshold- like changes due to the spin dependent interaction of cavity polaritons, are visualized. 
The observed dynamics cannot be reproduced within the conventional semi-classical model based 
on the Gross-Pitaevskii equations. To explain the observed phenomena, it is necessary to take into 
account the unpolarized exciton reservoir which brings on additional blueshift of bright excitons, 
equal in the and cr~ polarization components. This model explains the effect of polarization 
instability under both pulsed and continuous wave resonant excitation conditions, consistently with 
the spin ring pattern formation that has recently been observed under Gaussian shaped excitation. 



PACS numbers; 71.36.-|-c, 42.65.Pc, 42.55.Sa 

I. INTRODUCTION 

Cavity polaritons are bound light-matter states that 
appear due to the strong coupling of 2D excitons 
and photons in semiconductor microcavitiesjii^ The 
optically driven system of polaritons behaves like a 
weakly imperfect Bose gas, which results in a spec- 
tacular row of polariton collective phenomena such 
as multistability,— parametric scattering,—"— pattern 
format ion fSiiS self-organization effectSf^iii^ dynamical 
Bose-Einstein condensation,-- Polariton multistability 
attracts much attention as a potential candidate for very 
fast picosecond range optical switching on a micron size 
scale. ^^'^^ Moreover, the sharp transitions in intra-cavity 
field, stemming from the multistability, can proceed con- 
currently with inter-mode parametric scattering, which 
determines interesting ways of polariton self-organization 
under pumping near the inflection point of the lower po- 
lariton branch (the "magic angle" ) i ^'^-'-i-'-^ 

The polariton bistability has recently been studied in 
the "scalar" approximation neglecting the spin degrees of 
freedom. '^'^ Because of the mutual interaction of polari- 
tons, their energy effectively depends on the intra-cavity 
field magnitude. The repulsion between excitons leads to 
the blueshift of polariton energy. When the pump is it- 
self blue-detuned from the polariton resonance, a strong 
positive feedback loop between the resonance energy and 
the field amplitude is created in a certain range of sys- 
tem parameters, resulting in sharp jumps in the intra- 
cavity field and, hence, in the cavity transmission. Con- 
sequently, the transmission signal can exhibit a promi- 



nent hysteresis in the dependence on continuous-wave 
(cw) pump intensityj^i^ 

The response of the optically excited polariton system 
becomes more complicated with allowance made for ex- 
citon spin degrees of freedom (that corresponds to arbi- 
trary optical polarization of the intra-cavity field) . In the 
general case the system has up to four stable states un- 
der a given cw pump, whereas the actual state of the 
system is determined by the history of the excitation 
process4^ In critical points, where a number and/or 
stability of stationary solutions change, the system can 
exhibit sharp jumps in both the amplitude and polariza- 
tion of the intra-cavity field. Under a spatially inhomoge- 
neous (e.g., Gaussian shaped) cw excitation the system 
can also exhibit a nontrivial spatial distribution of po- 
larization of the luminescence signal, like the "spin ring" 
patterns,^'* due to the same underlying phenomena. 

The strong multistability effect predicted in Ref. 
has recently been observed experimentally)^ including the 
sharp jumps in the cavity transmission under a smooth 
variation of the pump polarization degree. The spin 
ring patterns under Gaussian shaped excitation have also 
been reportedi ^^i^^ All these experiments were carried 
out under the cw excitation. Thus, one of the still re- 
maining questions is the dynamical peculiarities of the 
transitions in a multistable system, which could only be 
traced using the time-resolved techniques. Particularly, 
the characteristic switching times between different sta- 
bility branches call for the experimental study, for they 
can have a crucial impact upon the practical implemen- 
tations. 
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Another question concerns the theoretical approach al- 
lowing to describe the observed multistability effect in 
microcavities. Traditionally, the multistability is consid- 
ered in terms of the self-acting classical fields correspond- 
ing to macro-occupied coherent polariton modes which 
appear under a coherent resonant excitation.— i^ii^ii^ Al- 
though such an approach is supposed to be sufficient to 
describe the bistability in a circularly polarized system,'^ 
it gives wrong predictions for the general case of ellipti- 
cally polarized excitationJ^ Most probably, the incoher- 
ent states of the exciton reservoir which are inevitably 
excited in optical experiments have a substantial impact 
on the decay rates and energies of polaritons. For in- 
stance, the nonlinear decay of polaritons with different 
circular polarizations was taken into consideration in or- 
der to explain the experimental data in Ref. jSt Further, 
the nonequilibrium transitions reported in Ref. [13 can 
only be reproduced in calculations taking into account 
the reservoir induced shifts of the polariton energy. 

In the present work we report the experimental study 
of nonequilibrium transitions in the multistable cavity 
polariton system. Unlike in recent Refs. [13, and [la, 
we have studied the optically polarized system under a 
pulsed nanosecond-long excitation, which allowed us to 
trace the time-resolved dynamics of the intra-cavity field. 
The shifts in polariton energy were reflected by temporal 
variations in the transmission energy spectrum. Thus, 
the employed technique is capable of visualizing the tem- 
poral correlations between the resonance energy and in- 
tensity of the intra-cavity field. 

We discuss in detail the time dependence of the trans- 
mission signal polarization, below as well as above the 
threshold, for several polarizations of the pump beam. 
The observed polarization behavior cannot be repro- 
duced within a semi-classical model based on the Gross- 
Pitaevskii equations considered in Refs. [3, [13; and [l^. 
The experimental results allowed us to develop a phe- 
nomenological model to describe a multistable polariton 
system with a proper regard to the exciton reservoir. The 
reservoir excitons shift the polariton energy that, in turn, 
influences the thresholds of nonequilibrium transitions. 
Both the linear and nonlinear mechanisms of exciton 
scattering into the reservoir are found to be significant 
within a sub-nanosecond time scale. While the nonlinear 
decay of cross-circularly polarized excitons leads to the 
levelling of the and a~ jump points (which has also 
been found in Refs. [^ and [l3 in the cw pump regime), 
the linear decay mechanism leads to the temporal delay 
of the jumps with respect to the pump intensity peak. 

Further, we apply the developed approach to simulate 
the parametric scattering under pumping at the magic 
angle. As we show, this model allows to explain self- 
consistently the temporal dependences of polarizations 
of both the driven mode and the scattering signal which 
appears at the polariton branch bottom. 

The paper is organized as follows. In Sec. [ll] the avail- 
able models of the multistability effect are considered and 
compared with each other in view of the recent experi- 



mental results of Refs. H, [H and M- Sec. ^ references 
the experimental setup. Sec. IIVI contains the experimen- 
tal results and compares them with the calculations per- 
formed in the framework of the suggested model, for the 
cases of pumping at normal incidence (Sec. IIV Ap and 
at the magic angle (Sec. IIVB|) . with a brief reference to 
some of the still unresolved issues fSec. lIV Ct . The results 
are summarized in Conclusion (Sec. |Vl) . 



II. THEORETICAL MODELS 

The multistability effects in microcavities were recently 
considered in the framework of Gross-Pitaevskii equa- 
tions written for the macro-occupied polariton modesi^ii^ 
These equations can also be written for the strongly cou- 
pled exciton and photon fields {V and £ respectively) in 
the cavity active layeri^^ 

i£+ = (c^c - He) £+ + aJ'+ + PV+, (1) 

i£- = (cjc - ilc) £- + aT- + PV- , (2) 

iV+ = {u^ + + V2\V-\'' - ij^)V+ + A£+, (3) 

iV- = (c^x + V2 |P+ 1' + Vi \V- 12 - ii^)V- + (4) 

(in the cr^ basis). Here, F stands for the incident electric 
field that is usually treated as a plane wave: F± oc e~*"p'; 
cjc.x and 7c, X are the eigenfrequencies and decay rates of 
the intra-cavity photon and exciton modes; a and /3 are 
the cavity response coefficients, and A is exciton polar- 
izability (so that 2^Aji equals Rabi splitting); V1.2 are 
the matrix elements of the interaction between excitons 
with same (Vi) and opposite (V2) circular polarizations. 

For simplicity, Eqs. ([THll) are written for only the 
driven polariton mode with zero quasi-momentum, al- 
though they can be easily generalized to the case of es- 
sentially many-mode systemJ^ In the latter case, the us- 
age of the exciton-photon basis (instead of the simplified 
polariton basis) brings the advantage of taking account 
of the TE/TM splitting effects as well as the angular de- 
pendence of the polariton interaction constants. As will 
be shown below, it is also convenient for taking account 
of the incoherent exciton reservoir. 

The multistability conditions are defined by the re- 
lation between exciton-exciton interaction constants, Vi 
and V2. According to the study of polariton spin quan- 
tum beats in a GaAs cavity,^° the cross-circularly po- 
larized excitons exhibit attraction (V2 < 0) that is 
weaker than the repulsive interaction of co-circular ex- 
citons {Vi > and 1^21 < ^i)- The fact that V2 is 
negative is also confirmed by the rotation of polariton 
polarization direction during the process of parametric 
scattering. In the general case the processes of exciton- 
exciton interaction are mediated by biexciton and dark 
exciton states,—"— hence the relation between the inter- 
action constants of polaritons can be dependent on the 
exciton-photon detuning. More precisely, one should 
also consider the saturation of excitonic transitions as an 
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additional nonlinear mechanism affecting the blueshift of 
polaritons under a sufficiently strong excitationi^ 

If Vi > and V2 < (as assumed in Refs. 0, [13, and 
[igI basing on the experimental estimations of Ref. [20| ). 
then the cr^ -components of the intra-cavity field are al- 
most uncoupled. If V2 < 0, then the threshold intensity 
W = l-F+P -|- 1-7^- P of linearly polarized pump is at least 
two times larger than that of circularly polarized pump: 

^(lin) > 2^(">^<=) 

Contrary to the above predictions, 



(lin) 



found to be slightly less than W^^^'^^ in the exper- 
iments using cw as well as time-resolved pumping 
techniques.— li^i^i^In order to explain it, two theoret- 
ical approaches have been proposed. The first one^i^^ 
presumes a repulsive interaction of polaritons with oppo- 
site pseudospins (cr^ polarizations) in the framework of 
Gross-Pitaevskii equations of type ©-(HI). The second 
approacbi2i2^ considers the additional blueshift resultant 

^(circ) 
thr 



thr 



has been 



from an exciton reservoir, which leads to W^^^^ < 
without the necessity for a repulsive interaction between 
cross-circularly polarized excitons. (A somewhat similar 
model has previously been proposed in Ref. 'so' to describe 
the interaction of a "spinless" polariton condensate with 
a reservoir.) 

In the present work we explore the second model 
accounting for the transitions of the optically driven 
excitons into the incoherent (reservoir) state in which 
the overall pseudospin is relaxed, so that the reservoir 
provides equal blueshifts for both polarization compo- 
nents of the coherent state.'"' ''■-^^ Those transitions are in- 
troduced phenomenologically as non-radiative decays of 
bright excitons accorded with the rate of occupation of 
the optically inactive reservoir. The model for the intra- 
cavity electric field £ and exciton polarization V coupled 
with the integral population of the reservoir Af can be 
written as 

i£+ = {lUc - ilc) £+ + aT+ + PV+ , (5) 
i£- - (wc - ilc) £- + olT- + PV- , (6) 
iV+ = [u^ + Vi\V+\^ + V2\V-\^ + {Vi + V2)N/2 - 

i{l. + l„ + V,\V-\^)]V++A£+, (7) 
iV^ = [u^ + Vi\V-\^ + V2\V+\^ + {Vi + V2)N/2 - 

*(7x + 7xr + K-im')]7'-+Af_, (8) 
N = -7rAA + 27,, (IT'+P + |P_|2) + 

m'P+?\V-\''. (9) 

This model is a generalization of Eqs. (IT])-(|1]). Here, 7xr 
is an additional decay rate of excitons that corresponds 
to the linear mechanism of light absorption, providing 
the term 2-f^,{\V+\'^ + \'P-\'^) in Eq. dH); V, stands for the 
rate of nonlinear interaction between V± that provides 
an additional occupation of the reservoir (4T/i.|'P+P|'P-P 
per unit time) due to the mixture of excitons with op- 
posite polarizations; 7r stands for the reservoir own de- 
cay ratejSi Microscopically, the nonlinear absorption of 



cross-circularly polarized excitons is resultant from the 
scattering of a pair of bright excitons with opposite spins 
[ Jz = —1 and Jz — +1) into dark excitons (J^ = —2 
and Jz — +2), which is closely related to biexciton cre- 
ation (see Ref. 22 for details). We do not consider the 
reverse transitions of incoherent (reservoir) excitons into 
the driven mode; this approximation is valid as long as 
the occupation of the reservoir states is small enough. 
Under resonant pumping at k = the energy mismatch 
of the pump with respect to free exciton makes the filling 
of reservoir quite ineffective, but it is still possible due to 
energy level broadening and non-zero temperature. The 
reservoir can accumulate excitation due to its compara- 
tively long lifetime (see Sec. IIV|) . 

Equations ©-([9]) can easily be generalized to the 
case of many-mode system, with exactly the same non- 
linear interaction terms (which are "local" in the real 
space), and a linear part being in the fc-space the same 
as in Ref. 1I61 . In the many-mode calculations repre- 
sented below in Figs. [21 |3l and [TT] [T3l the cavity dis- 
persion £^™'™(fc) and response coefficients aTE,TM(fc), 
/3TE,TM(fc) for TE and TM cavity modes are properly 
taken into account using the transfer matrix techniquci^ 

The two models, Eqs. ©-(gl) and ([5])-®, can give 
similar predictions for the cases of circular and linear 
polarizations of excitation. In particular, the key experi- 
mental effect waiting for explanation — that the threshold 
of linearly polarized pump is smaller (W^^"^ < W^^^'^'') — 
can be reproduced within the Eqs. dl])-© with V2 ^Vi- 
Nonetheless, the transitional dynamics under the inter- 
mediate pump polarizations < pp < 1 are completely 
different in the two systems, see an example in Fig.m 

If < and K- > (Fig. [Jh, solid line, and Fig. O), 
the cavity transmission polarization increases with pump 
in the sub-threshold area, due to the dominant blueshift 
of the leading polarization component (pp = -1-0.6). On 
reaching the threshold, the leading component jumps up 
to the high-energy state, which causes a step-like in- 
crease in reservoir population and the blueshift of cr~ 
mode. Consequently, the minor mode also enters the 
above-threshold area, and the degree of circular polar- 
ization (DCP) of transmission (ptr) decreases down to 
the values smaller than the pump DCP. The same ef- 
fect results in formation of the spin ring patterns under 
the Gaussian shaped excitationi^iiiii^ in a close vicinity 
of the threshold pump power (Fig. [5]). As a matter of 
fact, the ring shaped pattern of DCP indicates that both 
polarization components reach the high-energy state at 
the center of the pump spot where, therefore, the trans- 
mission DCP is minimal. The spatial distribution of ptr 
(Fig.[2lff-/) exhibits a well pronounced minimum at a; = 
already when the threshold is surpassed by 1.5 times 
(Fig. [2]e), in agreement with the experimental observa- 
tions (see Fig. 3 in Ref.flTt). This markedly differs from 
a prediction of the conventional theory based on Gross- 
Pitaevskii equationsiiiiii£ with \V2\ < Vi (see Fig. O 
which does not exhibit the rings even at W/Wthr ~ 5 for 
Pp = 0.6. The temporal dynamics for IV2I ^ Vi is shown 
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FIG. 1. (a) The calculated time dependence of the cav- 
ity transmission polarization under a slowly varying pump 
(transverse-dashed triangular line) with elliptical polariza- 
tion. The degree of circular polarization (DCP) of pump is 
Pp — 0.6 (thin solid line). Thick solid line shows the data cal- 
culated using Eqs. ©-(© with V2/V1 = -0.1, V^i-jrVi) = 
6meV~^, and 7xr = (explicitly as in Ref Dash- 
and-dot line represents the same dependence calculated using 
Eqs. ©-(HI) with V2/V1 = -0.1 (this ratio was used in Refs. 4 
and 1141'). Dash line corresponds to the case of V2/V1 = +1.1. 
The calculated data is time-averaged for each temporal point 
over 100 ps intervals, in order to eliminate the fast transi- 
tional effects which are not observable under the cw pump 
conditions, (b-d) The corresponding time dependences of the 
intra-cavity field intensity in polarization components. 



in Fig. [T]a (dash-and-dot line) and in Fig. [T]c. 

In the opposite case, when V2 > V\ and the system is 
"purely coherent" (Fig.[T]a, dash line, and Fig.[T]ff), which 

also corresponds to the right relation W^'' < wl^"'^\ 
the polarization dynamics is completely different. First, 
the transmission DCP decreases in the sub-threshold 
area, since an increase in the "leading" (t+ component 
provides a larger blueshift for cr" , so that the difference in 
cr^ intensities weakens with increasing pump. On reach- 
ing the threshold, ptr shows a minor step-like increase fol- 
lowed by the very weak changes with further increasing 
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FIG. 2. The calculated "spin ring" pattern formation in 
the framework of many-mode model Eqs. ((5])-(|9]) with K = 
0.6 • Vi, = -0.1 • Vi, and 7xr = (explicitly as in Ref. II). 
The pump DCP is +0.6 (as in Fig. [!}. Panels a-e show the 
spatial distribution of intensity of the intra-cavity field in cr^ 
and a~ polarization components (solid and dashed lines, re- 
spectively), for several values of pump power W; the pump 
distribution profile is shown by dotted lines. Wo roughly 
corresponds to the threshold excitation density. Panels A-j 
represent the corresponding spatial distribution of the intra- 
cavity field DCP. 
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FIG. 3. The calculated spatial distribution of the intra-cavity 
field for K = 7xr = and V2 = — 0.1-Vi. All other parameters 
coincide with those of Fig. [2] Unlike Fig. [2] almost no ring 
structure appears in the cavity transmission DCP even at 
W/Wo = 5. 



pump power. No prerequisites for the spin ring formation 
are satisfied in this case. 

Thus, the observed phenomena — the disparity of 
< M^t*hr'' together with the field 



thresholds, W^/^"^ 
polarization dynamics — cannot be reproduced within 
Eqs. using only two constants Vi_2 of exciton- 

exciton interaction even if a single polariton mode at 
k = is considered. The experiments revealing both 



the energy (blueshift) and polarization dynamics of the 
driven polariton mode under the nanosecond-long exci- 
tation pulses are described below in Sec. Hvl 
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III. EXPERIMENT 

The microcavity structure grown by a metal-organic 
vapor-phase epitaxy technique has top (bottom) Bragg 
reflectors composed of 17 (20) A/4 Alo.13Gao.87As/AlAs 
layers. The 3A/2 GaAs cavity contains six 10 nm thick 
Ino.06Gao.94As/GaAs quantum wells. The Rabi split- 
ting is about 6meV. A gradual variation of the active 
layer thickness along the sample provides a change in 
the photon mode energy Ec = fiujc and, accordingly, in 
the detuning A between the exciton Ex{k=0) and pho- 
ton Ec{k = 0) mode energies. Experiments are carried 
out in several regions of the same sample with A w 
(Sec. and A « -1.5 meV (Sec. HYbI) . 

The sample is placed into the optical cryostat with 
controlled temperature. To excite the cavity, we use a 
pulsed Ti:sapphire laser producing picosecond long pulses 
at a repetition rate of 5 kHz. Prior to coming into the 
cryostat, the pulses pass trough a long multi-mode opti- 
cal fiber and then through a monochromator. After the 
fiber, the pulses have a duration of about 1 ns, and af- 
ter the monochromator they have spectral full width at 
half maximum (FWHM) of 0.7 meV. The excitation is 
performed either along the normal to the cavity plane 
(fcp = 0) or at the magic angle (fcp = 1.8 ^m~^) slightly 
above the lower polariton branch resonance, with var- 
ious pump polarizations (circular, linear, or elliptical). 
The main axis of pump polarization is directed along the 
(110) axis of the structure. The pump beam is focused 
onto the spot with a diameter of 100 /im. The kinet- 
ics of the cavity transmission signal /tr(i) is detected at 
T « 6 K by the streak camera with spectral, angular, and 
time resolution of 0.28 meV, 0.5°, and 70 ps, respectively, 
in various polarizations. The temporal dependences of 
the signal are averaged over ~ 10^ pulses, for we have to 
collect the signal in order to eliminate noises. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 

A. Pumping at normal incidence (fcp = 0) 

The information on the intra-cavity field kinetics is ob- 
tained from the cavity transmission measurements J^i^^ 
The active region of the cavity is separated from a de- 
tector by the Bragg mirror that does not introduce any 
nonlinearity and/or spectral selectivity. Thus, the in- 
tensity of the transmission signal, Iti{(^,t), is propor- 
tional to the squared magnitude of the intra-cavity elec- 
tric field |£(kp,w,i)p, and the first momentum E — 
h J ujItT:{uj, t) duj/ J ItT:{uj^ t) doj at k = kp refiects the time 
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FIG. 4. Time and energy dependences of the excitation pulse 
(a, e) and cavity transmission intensity Itr = Itr{E,t) for dif- 
ferent pump powers W {b-d and f-h). The pump is polarized 
circularly. The normalized temporal shape of the excitation 
is shown at the left side (white solid line); the weighted mean 
value of the pump energy is shown by the vertical line. Left 
and right panels represent the experimental and calculated 
data, respectively. The color scale is logarithmic. 



dependence of the average (effective) energy of excited 
mode E{t). 

Figures H^^-d show the measured time dependences of 
transmission signal Itrit) under the excitation with cir- 
cularly polarized (cr"*") light. The spectra are recorded in 
the same cr"*" polarization; the fraction of signal in the 
opposite polarization {a~) does not exceed a few per- 
cent. At low peak excitation density W = maxj Itr{t) ~ 
9.5kW/cm^ the maximum of the signal is found at the 
low energy side of the excitation, near the lower polari- 
ton eigenfrequency. Although the blueshift in polariton 
energy remains relatively small, the time dependence of 
the signal differs markedly from that of the pump pulse. 
This indicates the onset of nonlinear processes at the red 
edge of pump spectrum. In particular, the maximum of 
the signal is delayed for a few hundred picoseconds with 
respect to the pump peak. Nonlinearities in both the 
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Time (ns) 

FIG. 5. Time dependences of tlie cavity transmission inten- 
sity, in circular and linear polarization components (cr* and 
TTx.y, resp.) for different pump polarizations pp (circular, el- 
liptical, and linear polarizations). The top panels show the 
time shape of the exciting pulse. Left and right panels rep- 
resent the experimental and calculated data, respectively. In 
the experiment, W = 28.5kW/cm^. 



magnitude and spectral position of the signal increase 
with pump power: as the peak energy of the signal shifts 
up to the pump energy Ep, (i) the signal intensity shows 
the superlinear (threshold like) increase and (ii) the time 
when the signal reaches its maximum value shifts back- 
ward to the onset of excitation pulse. 

The observed correlation between the energy blueshift 
and intensity of the signal is expected in the framework 
of the multistability model^ Eqs. however this 

model fails to explain the process of instability develop- 
ment. Figures ^-h show that the observed time and 
spectral dependences can be qualitatively reproduced in 
simulations using Eqs. (IS|)-(jni) with appropriate time and 
energy shapes of the excitation (see Appendix) . The cho- 
sen parameters of microcavity coincide with those used in 
the experiment; the constants of interaction with reser- 
voir are considered below. 

The temporal dependence of transmission under the 
excitation with linearly, elliptically and circularly po- 
larized pulses (DCPs 0, 0.2, 0.6, and 1) at = 
28.5kW/cm^ is shown in Fig. [SJ The temporal delay 
of the dominating cr"*" and tt^; polarization components of 
the signal with respect to the pump peak is seen to be 
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FIG. 6. Time dependences of degrees of circular {b, e) and 
linear (c, /) polarizations of the transmitted pulses, for sev- 
eral values of the pump DCP {pp = 0.2,0.4,0.6). Left and 
right panels represent experimental and calculated data, re- 
spectively. In the experiment, W = 28.5kW/cm^. 




Energy (meV) 

FIG. 7. The spectrum of intra-cavity filed in cr* compo- 
nents in the range of polarization instability (experiment). 
Weighted mean values of the energy E± are indicated by ar- 
rows. 



in the range of 150- 200 ps, which is much larger than 
the lifetime of polaritons with k = 0. The "anomalous" 
increase in the intra-cavity field observed in the range 
of decreasing pump power (after the pump peak is al- 
ready passed) can be explained by the blueshift in effec- 
tive polariton energy. In turn, this blueshift could only 
be provided by the long-lived polariton states accumu- 
lated during the pulse. Such an effect corresponds to the 
"linear" process of polariton scattering into the exciton 
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reservoir, which is determined by the ratio between 7r 
and 7xr in Eqs. ([T])-®- The experimental results are 
qualitatively reproduced in numerical simulations with 
= -0.1, V^lVx = 7 • 10-3, hjr ^ 2 - 10-3 meV, 
and 7xr — 3 • 10"^ meV. These parameters are chosen to 
meet the following conditions: (i) the overall occupation 
of the reservoir is comparable to that of the driven (op- 



(circ) 



thr 



tically active) polariton mode, so that W^^^' 
and (ii) the temporal peak of the signal can be delayed 
by hundreds of picoseconds with respect to the excitation 
peak even for pp — ±1. Note the long decay times of the 
reservoir states, ^ 300 ps. So long times are characteris- 
tic of excitons localized due to fluctuations of quantum 
well potential and/or free excitons with large lateral wave 
numbers. 

Figures and [S]e reveal the transmission signal to 
be retaining both circular and linear pump polariza- 
tions; in these cases only a weak depolarization of the 
transmitted light is observed. Retaining of the circu- 
lar polarization in agreement with the angular momen- 
tum conservation law is expected. Measurements per- 
formed under linearly polarized pump in a wide range 
oi W — 9.5 — 28.5 kW/cm'^ (not shown) have revealed 
the magnitude of transmitted light depolarization to be 
within the limits of 12-17%, almost independently of the 
pump power. On the other hand, under the elliptically 
polarized excitation (0 < pp < 1) the signal DCP does 
depend on the pump intensity due to the anisotropy in 
polariton-polariton interaction)^ Accordingly, the ratio 
between both circularly and linearly polarized compo- 
nents of the signal intensity (cr^ and TTx,y, respectively) 
varies with time (Fig. [S]c,(i). The corresponding tem- 
poral dependences of the degrees of circular and linear 
polarizations of the signal are shown in Fig. [6] From 
comparison of Figs. [S] and [HI it is seen that the reaching 
of maximum signal intensity at t ~ 0.2 ns is accompanied 
by a decrease in the signal DCP, in much the same way 
as the cw driven system behaves (Fig.[T]), though the lat- 
ter system exhibits the expectedly sharper transitions. 
On the other hand, our experimental technique allows 
us to visualize the energies of intra-cavity field in addi- 
tion to their intensities. The information on the energies 
is extracted from the transmission spectra as shown in 
Fig. [7] since Itr{huj) cx \£{hij)\'^. The energies of intra- 
cavity field in ct+ and polarizations demonstrate a 
well pronounced blueshift with increasing polariton den- 
sity and a well pronounced splitting. In particular it is 
seen in Fig. [7] that in the range of maximum polarization 
the component of the field is markedly (by 0.3 meV) 
blue-shifted with respect to . 

The basis in which shown degree of linear polarization 
(DLP) was measured is related to the pump polarization 
axis. We observe no linear polarization in the basis ro- 
tated by 45° with respect to the former one. Thus, the 
sum of squared Stokes parameters per signal intensity 
is less than 1, which means that the signal is partially 
depolarized; this is because the signal polarization com- 
ponents are measured as averaged quantities. 
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FIG. 8. Time dependences of the transmission DCPs, for 
several values of the pump power W, under excitation with 
elliptically polarized pulses with pp — 0.4. Grey dotted line 
corresponds to the normalized temporal shape of excitation. 
Left and right panels represent experimental and calculated 
data, respectively. 



In order to clarify the physics underlying the signal po- 
larization transitions, let us consider the series of mea- 
surements performed under elliptically polarized pulses 
with Pp = 0.4 and various pump intensities. Fig. [S] rep- 
resents measured (left side) and calculated (right side) 
time dependences of the transmission signal in tr^ polar- 
ization components. From these series, the pump power 
dependence of the temporal delay between the signal and 
pump peaks is evidenced. The certain conditions for a 
switch in signal polarization depend on both the pump 
intensity and its temporal shape. In case of the very low 
pump densities, the intra-cavity field is insufficient for 
multistability to reveal; on the other hand, if the pump 
density is increasing quickly enough then a jump into a 
high-energy state happens with no delay related to the 
long-lived reservoir, since most of the intra-cavity field 
is concentrated in the driven (optically active) polariton 
state. Except the case of maximum W, the dynamics 
shown in Fig. |5] lies in the intermediate range of pump 
powers which all are insufficient to pull a polariton state 
by itself up to the threshold magnitude. Instead of that, 
the non-equilibrium transitions proceed through the me- 
diation of excitation accumulated by the reservoir. In 
such a case the development of instability in the leading 
polarization component (cr+) always precedes the reser- 
voir induced switch in signal polarization. 

Figure [HI which also shows a series for two peak pump 
intensities at fixed pp = 0.4, allows us to visualize the 



8 



1.0 



c 
o 

.2 Ph 



0.5 



0.0 



0.6 



0.4 




a 

I Q 0.6 



0.4 

-1.6 

-2.0 
-2.4 
-2.8 



13.3kW/cm2 
19.0 kW/cm^ 
pump 




(c) 







— 25 arb.units 
32 arb.units 

— pump 



(f) 





0.0 0.5 1.0 0.0 
Time (ns) 



0.5 



1.0 



FIG. 9. Time dependences of the pump pulse intensity (a, e), 
degrees of circular (6, /) and linear (c, g) polarizations (DCP 
and DLP) of the transmitted pulse, and the weighted mean 
values E± of the transmission energy in (thick lines) and 
a~ (thin lines) polarization components {d, h), for several 
values of the pump power W. The pump DCP is 0.4. The 
mean value of pump energy, which exhibits minor variations 
in a close vicinity oi E — = —2 meV, is shown by solid 
thin lines in panels {d, h). Left and right panels represent 
experimental and calculated data, respectively. 



temporal correlations between the energy (in both cr* 
polarization components) and polarization degree of the 
transmission signal. The overall process of the instability 
development is accompanied by blueshifts in cr^ energies. 
At the first stage (which starts the earlier, the higher is 
pump density W), the leading polarization component 
(cr"*") enters the region of positive feedback between its 
effective resonance energy and the cr'^ field ampli- 
tude, so the difference in ct^ energies grows along with 
the DCP of transmission signal (Fig. [S?;, d). Accord- 
ingly, a well pronounced decrease in the DLP of transmis- 
sion is observed at the same time interval (Fig. [S]c). At 
approximately the half way through the energy growth 
(see panel rf), the integral occupation of reservoir be- 
comes sufficient to provide the shift in the minor com- 
ponent (cr~), which leads to decrease in both the en- 
ergy difference — E- and DCP, and to the corre- 
sponding increase in DLP. All these effects are quali- 
tatively reproduced in calculations shown in Fig. [S]e-/i. 
Thus, the observed dynamics confirm the effect of reser- 
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FIG. 10. The effective delay of the a components of the 
signal with respect to the pump vs. pump power (as defined 
by Eq. ([10])), for two pump polarizations, pp = 0.2 (panels a, 
c) and pp — 0.4 (panels b, d). Left and right sides of the figure 
represent experimental and calculated data, respectively. 



voir, though the energy broadening makes the transi- 
tions much smoother compared to the cw driven system 
(Fig.©. 

Figure \W\ shows the time delay of the signal with re- 
spect to the excitation pulse in a wider range of pump 
powers W; it summarizes the consideration of the thresh- 
old effects in the cavity transmission. The value of the 
delay At is determined as the difference between the mo- 
ments the signal and the pump reach their half-maxima 
for the first time: 



AU 



/±(t)-^imaxt/±(t) 



/p(t)— i-^ maxt /p(t) 



(10) 



At small W the response is linear, and the delays At are 
expectedly small. At threshold, the superlinear growth 
of the signal is observed with At of about 150 to 200 
picoseconds, i. e. at the back front of excitation. Thus, 
this transition is induced by the reservoir as explained 
above. With further increasing W, the optically active 
fraction of excitons grows, which lowers the system's in- 
tertness, and the delay of response decreases. Note, that 
the minor (cr") polarization component exhibits the su- 
perlinear growth at the same peak power W as cr^ , and 
the relative delay of and a~ components, Ai+ — Ai_, 
remains nearly constant above the threshold. This is an 
extra evidence of the fact that the jump in the minor 
component is provided by the reservoir even at high W. 
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B. Pumping at the magic angle (fcp — 1.8 ^) 

Let us now consider the polarization properties of 
the signal of parametric scattering — generally referred to 
as cavity "optic parametric oscillator" (OPO) signaP — 
which appears at the branch bottom (fc ~ 0) if pumping 
is at the magic angle {kp « 1.8 /im~^). There are two 
points to discuss here. First, in the general case of el- 
liptically polarized excitation the signal is expected to 
reflect variations in the DCP of the driven mode that, in 
turn, are shown to be strongly affected by the excitonic 
reservoir. Second, if the pump is polarized linearly, the 
experiments show the 90° rotation of the signal polar- 
ization axisj^ thus, we have to check if Eqs. ([S])-® can 
reproduce this effect in spite of the reservoir induced in- 
teraction between the cross-circularly polarized excitons. 

The experiments are carried out at negative detuning 
between the photon and exciton (A « — l.SmeV) in or- 
der to have the pump frequency close to that employed 
for the normal incidence pumping considered in Sec IIV Al 
The cavity and excitation parameters used in the simu- 
lations coincide with the experimental ones, whereas the 
nonlinear interaction constants as well as reservoir char- 
acteristics are the same as used in the previous Section. 

Figure [TT] represents measured and calculated dynam- 
ics of intra-cavity field at fc = fcp (which corresponds to 
the cavity transmission) and fc (OPO signal) un- 
der elliptically polarized excitation with DCP = 0.5 and 
different pump powers W. When W exceeds the criti- 
cal value (in the experiment Wthi ~ lOkW/cm^), a well 
pronounced superlinear growth of the cavity transmis- 
sion (Fig. [Tit )) is observed at t = 0.1 -0.4 ns. Near the 
threshold, a doubling of the pump power W leads to an 
order-of-magnitude growth of the OPO signal intensity 
(Fig. [TT]s). The observed nonlinearities indicate the on- 
set of parametric instability, being, however, by far not 
as sharp as under cw excitationj^i 

The dynamics of the driven mode is similar to that 
under the pumping at normal incidence. The growth 
of the field amplitude, which continues with increasing 
pump density fFig. \TVb). expectedly leads to the growth 
of transmission DCP fFig.lTTb) due to the dominantly re- 
pulsive exciton-exciton interaction. Att = 0.1 -0.4 ns the 
transmission starts to behave strongly nonlinearly with 
respect to the driving field, and the peak of transmission 
becomes significantly delayed with respect to the peak of 
pump pulse. The latter indicates that the transition is 
assisted by the reservoir induced blueshift. Accordingly, 
the jump in transmission DCP is immediately followed 
by dropping back as soon as the system enters the above- 
threshold region where the energies of components of 
intra-cavity field (Fig.lTTIj) become partially levelled due 
to the reservoir filling. 

The domination of the leading polarization component 
(ct+) is strongly enhanced in the OPO signal (Fig. [TTf) 
whose DCP reaches 0.85 at t « 0.2 ns. The subsequent 
lowering of DCP of the driven mode forces the OPO sig- 
nal, too, to partially loss its polarization. However, the 
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FIG. 11. Time dependences of the pump pulse intensity (a, 
g), transmission intensity (6, h), transmission DCP (c, i), 
weighted mean values of the transmission energy E± {d, j) in 
two polarization components and a~ (thick and thin lines, 
resp.), OPO signal intensity (e, k) and DCP (/, I), for several 
values of the pump power W . Straight lines in panels c, i 
indicate the mean value of pump energy. Left and right panels 
represent experimental and calculated data, respectively. 



signal restores a high DCP during the further evolution 
accompanied by the decrease in transmission intensity. 
The calculations performed on the basis of Eqs. ([S])-®, 
which are presented in the right panels of Fig. Ill[ repro- 
duce the observed threshold-like behavior of intra-cavity 
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FIG. 12. Time dependences of tlie OPO signal intensity (a, 
c) and degree of linear polarization {b, d) under the pump 
^{t) polarized linearly along the TM cavity mode \\ £^tm), 
for some different values of the pump power W. Left and 
right panels represent the experimental and calculated data, 
respectively. The temporal shape of excitation is indicated in 
panels a, c as a guide for eye. 



field along with its polarization properties. However, a 
quantitative agreement is not achieved, partially due to 
a complicated spectral shape of the pump pulses (see 
Fig. lU and due to intrinsic limitations of the suggested 
theoretical model (see discussion in Sec lIVCj) . 

Let us turn to the case of exactly linear polarization of 
the excitation at fcp > (Figs. [T^ and [T^ . In this case 
the OPO signal depends on whether the pump excites a 
pure cavity state — one of TE and TM modes which are 
the eigenstates of 2D photon in the empty cavity — or a 
mixture of both the TE and TM components. Fig. [TH 
represents the dynamics of OPO signal under the pump 
polarized along the TM cavity mode. In the range of high 
intensity of the signal, its degree of linear polarization 
(DLP) reaches 90%. The polarization axis is 90° rotated 
with respect to the pump (what corresponds to negative 
DLP values), which is well reproduced in the calculations. 

On the other hand, if the pump is polarized in between 
the TE and TM directions {f \\ [^te+^tm], Fig.[I31), the 
signal DLP does not exceed 30% in the range of a strong 
signal (Fig. [T51e). The shown DLP values were measured 
(and calculated) in the basis of directions parallel and or- 
thogonal to pump polarization; there is no marked DLP 
in the TE-TM basis as well. This effect is explained by 
the lifted degeneracy (TE-TM splitting) of cavity modes 
which leads to the variation with time of the tr^ phase 
shift. Accordingly, it causes the misbalance of and a~ 
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FIG. 13. Time dependences of the transmission DCP (a, /), 
transmission DLP {b, g), signal intensity (c, h), signal DCP 
{d, i), and signal DLP (e, j) under linearly polarized pump 
wave T II {£te + ^tm). Left and right panels represent ex- 
perimental and calculated data, respectively. The temporal 
shapes of excitation are indicated in panels a, c as a guide 
for eye. In calculation, the TE/TM splitting of the driven 
polariton mode equals to 0.06 meV. 



components of intra-cavity field, i. e. a reduced DLP and 
non-zero DCP of the driven mode (Fig. [T5b . 6), which 
is seen in both measured and calculated data. Due to 
permanent modification of the direction of signal polar- 
ization, the eventual OPO signal appears to be signifi- 
cantly depolarized (the sum of squared DCP and DLP is 
markedly less than 1, Fig. [T5H. e). 

Thus, the dynamics of the OPO signal confirms the ef- 
fect of reservoir. Under elliptic excitation, the cr^ polar- 
ization components tend to be levelled by the reservoir. 
By itself, this effect could be assigned to a repulsion be- 
tween cross-circularly polarized excitons (V2 > 0). On 
the other hand, the axis of signal polarization is found in- 
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verted with respect to the pump, which restricts the sign 
of V2 to be negative irrespectively of the reservoir induced 
blueshifts. Since we have reproduced these effects within 
Eqs. ([S])-®, it is proven to provide a self-consistent de- 
scription of the intra-cavity field polarization dynamics. 



C. Shortcomings 

We have shown that the semi-phenomenological 
Eqs. ©-(in]) explain the observed polarization dynamics 
in a good qualitative agreement with the experimental 
data. Let us now consider the limitations of the sug- 
gested approach. Its most important shortcoming is the 
absence of a feedback between the driven polariton modes 
and the exciton reservoir, as long as the latter is intro- 
duced using the integral occupation number Af, neglect- 
ing the actual microscopic distribution of reservoir states. 
As a result, Eqs. ([1])-® do not describe, e. g. the estab- 
lishing of equilibrium between coherent and incoherent 
exciton "phases" which would normally take place in a 
real system; this process implies (i) the reverse transi- 
tions from the reservoir into the macro-occupied mode 
(similar to Bose-Einstein condensation) and (ii) depen- 
dence of the scattering rates (7xr and K-) on the polariton 
density. These shortcomings lead to several visible dis- 
crepancies between the measured and calculated data. In 
particular, the superlinear drop in transmission intensity 
observed at the back front of exciting pulses, which is ev- 
ident from comparison of left and right panels of Fig. |8l is 
not reproduced by the present model. Note as well that 
the nonlinear scattering rate should depend on the 
energy of biexciton resonance which mediates the transi- 
tions of optically driven polaritons into the reservoir . ^^1^^ 
Accordingly, the actual relation between the multistabil- 
ity thresholds W^^"'^'' and VF/h™^ would depend on the 
pump frequency as well as on the frequency detuning be- 
tween photon and exciton modes. 

As a result, we expect Eqs. ([S])-® to be qualitatively 
valid in the case of resonant excitation of the optically ac- 
tive cavity states. On the other hand, in the cases of com- 
parably large pump wave numbers and/or non-resonant 
excitation conditions, when a much larger contribution 
of reservoir excitons is expected, the model is likely to 
become unsatisfactory. 



V. CONCLUSION 

In the present work we have studied the non- 
equilibrium transitions in a multistable system of cav- 
ity polaritons under resonant nanosecond-long excita- 
tion. Using the spectrally broadened pulses, we have vi- 
sualized the temporal correlations between the effective 
resonant energy, intensity, and optical polarization of the 
intra-cavity field which all undergo the strong changes 
on reaching the threshold pump power. In the vicinity of 



the threshold, the dynamics of such system is strongly af- 
fected by the long-lived exciton reservoir (excited due to 
polariton scattering) which influences both the character- 
istic times of instability development and output signal 
polarization. 

The temporal behavior of the intra-cavity field is found 
to be not described in the conventional model based on 
the Gross-Pitaevskii equations written for purely coher- 
ent macro-occupied polariton states. Most importantly, 
the observed phenomena cannot be explained even quali- 
tatively within a model with only the two exciton-exciton 
interaction constants (Vi_2) allowed for. To explain the 
experiments, we have proposed the model for the macro- 
occupied polariton states coupled with an exciton reser- 
voir. In spite of some limitations, this model provides 
a self-consistent description of the observed intra-cavity 
field dynamics under both pulse and continuous wave ex- 
citation conditions. 

Authors thank M. S. Skolnick for rendered sam- 
ples, D. N. Krizhanovskii for fruitful discussions, and 
A.V. Larionov for assistance in the experiment. This 
work was supported by Russian Foundation for Basic Re- 
search and the Russian Academy of Sciences. 



Appendix: Details on numeric modeling 

In attempt to reproduce the experimental pump spec- 
trum, which is characterized by FWHM of about 0.8 meV 
and, hence, has a lowered coherence, the modeled pump 
is 'broadened' by means of the partial randomization of 
phases of different spectral harmonics. Namely, the pump 
amplitude is defined as follows: 

J-±(t) = J-f (i) exp(-^..p^ - z^±) ; (A.l) 

7^(t) = 5Iexp -^-*w„i-^^X™ . (A.2) 

n— 1 \ m— 1 / 

Here, 

• J-^^^ (t) are real- valued amplitudes of the polar- 
ization components; 

• — defines the direction of the polarization 
axis; 

• ujp is the central pump frequency; 

. w„ = r!(f - i), n^l,2,...N; N ^ 1000; n = 
2meV; 

m a = FWHM/V8 log 2 defines the width of the 
pump spectrum, FWHM = 0.8 me V; 

• {Xm I 'Ti = 1,2,...,A^} is a set of random uni- 
formly distributed numbers within the interval 
[0;1). 
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The spectral FWHM of the pump defined such a way 
approximately coincides with experimental one, and its 
coherence time is about 5 picoseconds. Under these con- 
ditions the system is inevitably stochastic. Hence, for 
each set of parameters Eqs. ([I])-® are solved a number 
of times (20-100) with different {x} until the averaged 
evolution of the signal intensity becomes independent on 



the number of realizations. This corresponds to the col- 
lection of numerous signals, which allows us to reduce 
noises. The degrees of polarization are then calculated 
using the averaged intensities. As a result, the transmis- 
sion can be partially depolarized (under elliptically po- 
larized excitation) due to small random differences in the 
individual pulses in much the same way as the measured 
signal. 
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For the sake of simplicity, only the nonlinear mechanism 
of reservoir population (7xr = and Vr 7^ 0) has been 
considered in Ref. [itI . In such a case, the reservoir in- 
duced blueshifts are determined mainly by the ratio be- 
tween reservoir's decay rate and the coupling constant V^. 
For instance, the calculations performed with different 7r 
but constant ji/V,: and a sufficiently slow variation of the 
pump power lead to exactly same results for the cases of 
circular and linear pump polarizations. In this connection, 
the experiments carried out with time resolution using the 
nanosecond long pump pulses (described in Sec. lIVp are ca- 
pable of clarifying the characteristic lifetimes of the reser- 



voir states, as well as the relation between linear and non- 
linear mechanisms of polaritons' scattering into the inco- 
herent exciton reservoir. For instance, the former mecha- 
nism is expected to play a role irrespectively of pump po- 
larization, affecting the dynamics of polariton energy and 
the time moments when the nonequilibrium transitions oc- 
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